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Abstract 

High Pj measurements of hard hadrons or jets at RHIC and LHC appear contradictory and in some cases counter- 
intuitive, but upon closer investigation they represent a coherent picture of jet-medium interaction physics which 
can be established with systematic comparisons of models against a large body of data. This picture is consistent 
with a perturbative QCD mechanism and does not require exotic assumptions. This overview outlines how several 
key measurements each partially constrain shower-medium interaction physics and how from the sum of those the 
outlines of the mechanism of jet quenching can be deduced. Most current jet results from LHC can be naturally 
understood in this picture. A short summary of what can be established about the nature of parton-medium interaction 
with current data is given in the end. 
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1. Introduction 

High transverse momentum (Pj) processes in Quantum Chromodynamics (QCD) have been suggested as a probe 
to investigate the physics of ultrarelativistic heavy-ion (A-A) collisions more than a decade ago and have become 
a cornerstone of the experimental heavy-ion program both at RHIC and LHC. Yet, compared with for instance the 
modelling of the bulk medium evolution using hydrodynamics where an era of precision fits to differential data has 
begun (see e.g. [ 1 1), theoretical progress in the physics of hard probes is slow. 

Partially, this can be attributed to the structure of the problem: Hard probes are at the same time sensitive to the mi- 
croscopical degrees of freedom of the medium and the macroscopical medium density distribution and evolution. Thus 
high Pj observables have the capability to potentially do both tomography and transport coefficient measurements in 
the medium, i.e. to constrain its microscopical and macroscopical dynamics. While the microscopical physics of the 
medium and hence the basic structure of the parton-medium interaction is currently in detail unknown, the density 
evolution (in terms of fluid-dynamical modelling) is constrained by bulk data, but not unambiguously determined. 
As a result, the description of any single observable in terms of a parton-medium interaction scenario and a medium 
evolution model is usually far from unique and does not allow a firm conclusion with regard to either question. 

However, systematic multi-observable studies of a matrix of parton-medium interaction and medium evolution 
models can reduce the allowed model space significantly, resulting in a class of models which shows good consistency 
across RHIC and LHC observables. 
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2. A classification of models 

2.1. The vacuum baseline 

In vacuum, at the core of the computation of hard processes in perturbative QCD (pQCD) is the factorization of 
the hard, short-distance interaction part and the long-distance physics of initial state parton distributions and final state 
parton showers. In a medium, QCD factorization is not strictly proven, but remains a plausible working assumption. 
A sufficiently large separation of hard and thermal momentum scale then ensures that the hard process itself receives 
no medium corrections and can serve as a 'standard-candle', producing high pj partons at a known and calculable 
rate. 

The partons emerging from the hard process are typically characterized by a large virtuality scale Q. If one takes 
the uncertainty relation to estimate the timescale at which intermediate shower states at lower virtuality can be formed, 
one finds t ~ E/Q 2 with E the virtual parton energy. Comparing this scale with typical medium lifetimes (9(10) fm, 
one finds that for typical RHIC and LHC kinematics a significant part of the partonic shower evolution (i.e. with 
Q > Aqcd,iiiii) takes place in the medium, with the shower being gradually boosted out of the medium at top LHC 
energies and on the other hand hadrons with large mass m/ r being produced in the medium. However, over a large 
kinematical range and for sufficiently light hadrons it can reasonably be assumed that the medium modification largely 
concerns the evolution of a parton shower with subsequent hadronization as in vacuum. 

As a baseline example for in-medium shower modelling, consider a virtuality-ordered parton shower such as 
implemented in the PYSHOW algorithm J2]. This is treated as an interated series of splittings of a parent into two 
daughter partons a — > be where Ef, - zE a and E c = (1 - z)E a and where the virtuality of partons in terms of 
t = In Q 2 /A 2 JCD decreases in each branching. The splitting probability is given by 



where P a ^bc(z) are perturbatively calculable objects, the so-called splitting kernels. The splittings cease when 
the parton virtuality reaches a non-perturbative lower scale 0(1) GeV at which point a hadronization model needs to 
be used. In this way, the original high virtuality of the shower initiator is converted into the transverse momentum 
distribution of the final shower partons (and the initial virtuality sets a kinematic bound for observables like jet shapes). 

An important thing to note is that the splitting kernels P a ^bc(z) are scale invariant, i.e. they do not depend on an 
absolute momentum scale. As a result, the fragmentation functions generated by a series of splittings is self-similar 
(the momentum distribution of hadrons in a subjet looks almost the same as the distribution in the whole jet when 
plotted as a function of z) and up to logarithmic corrections has no strong energy dependence. 

Estimating the timescale associated with the initial, hard branchings which largely determine jet shape and subjet 
structure, one finds r ~ 0.01 fm, i.e. well before any medium can be formed. This gives rise to the expectation that 
subject structure and jet shape at high Pj should not be very sensitive to medium modifications. 

2.2. In-medium showers 

There are two main processes by which this pattern can be modified in a medium: 1) direct loss of energy (by 
elastic collisions, drag, . . . ) into the medium as parametrized by a transport coefficient e and 2) enhanced soft gluon 
radiation by medium-induced virtuality as parametrized by a transport coefficient q. The first mechanism is in general 
incoherent and the energy is lost into the medium, the second mechanism involves a de-coherence time for the radiated 
gluon and transfers energy from the leading parton into an increased number of subleading shower partons. Note that 
this distinction can be made on average on the basis of transport coefficients, but not for a single interaction graph with 
the medium which may involve both direct and radiative energy loss. A recently suggested additional mechanism is 
the modification of color flow by the medium, leading to a modification of the hadronic shower evolution (3) . 

Models may now be classified based on what part of the dynamics they include, how they treat the medium and 
how they solve the shower evolution equations. The most complete approach are models which simulate the whole 
in-medium shower evolution such as JEWEL |4|, YaJEM |5 6], Q-PYTHIA |7| or resummed higher twist (HT) |8|. 
However, if one is only interested in observables which probe the leading shower fragments, such as single inclusive 
high Pj hadron suppression, which probe dominantly showers in which a large fraction of the energy flows through 
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a single parton, then leading parton energy loss modelling in which the medium evolution is cast into an energy loss 
probability P(AE) before vacuum fragmentation is applied is sufficient. Models of this type include the ASW |9|, 
AMY iflOl [TTI or WHDG lfl2l frameworks. Finally, in hybrid approaches a vacuum shower is computed down to 
some scale, then the produced partons are put on-shell and each propagated through the medium using a leading- 
parton energy loss model. An example for hybrid modelling is MARTINI lfT3"ll . 

Of considerable interest is how the medium is included into the modelling process. A number of models come 
with an explicit treatment of medium partons and their interaction with shower partons, typically as a thermal gas 
of quark and gluon quasiparticles J4] [TUl [TT| [T3l . Such a description requires that the medium is, at least to some 
degree, perturbatively tractable, which is not in line with the basic assumption underlying almost ideal hydrodynamical 
modelling which assumes that the medium is strongly coupled and the mean free path essentially vanishes. Usually 
an explicit treatment of interactions with medium partons breaks the scale invariance of the fragmentation function 
around the momentum scale of the medium partons and below, while it largely remains intact above. In addition, the 
medium provides additional kinematic phase space for transverse shower broadening. 

A different approximation is to cast all medium effects into a modification of the splitting probability, in essence 
replacing P a ^bc(z) - > P' a ^bS z ^ t ^ le detailed form of P' motivated by a specific physics scenario as done e.g. in 
Q-PYTHIA Q. Such a prescription explicitly conserves energy and momentum inside the shower, i.e. there is no 
additional phase space for transverse broadening. It also perserves scale invariance of the fragmentation function, but 
it alters its shape to a different functional form. 

Finally, the medium can be (without any explicit reference to its degrees of freedom) appear via transport coeffi- 
cients such as q, e and alter the kinematics of propagating partons as implemented in YaJEM |5][6]. This again breaks 
scale invariance of the fragmentation function around the medium momentum scale and below and leads to additional 
phase space for transverse broadening. 

There are two main strategies to solve the shower evolution equations — analytical and Mote-Carlo (MC) tech- 
niques. In general, analytical techniques get exact treatment of quantum interference effects while they often rely 
on kinematic approximations such as eikonal propagation or infinite parent parton energy which violate energy- 
momentum conservation. In contrast, MC frameworks usually have exact kinematics, but have to resort to phe- 
nomenological prescriptions for interference effects. 





medium modifies kinematics probabilistic picture explicit medium partons 


MC shower 
analytical shower 
MC hybrid 
analytical energy loss 


YaJEM OSl Q-PYTHIA Q, BW Ha JEWEL @ 

— HT resummed [8] — 

— PYQUEN El MARTINI 03 

— ASW GO, HT ED AMY IfTOlfTTl. WHDG 113 



Table 1 : A summary table of a rough classification of several well-known jet quenching models. 



These different approaches to solving the in-medium shower evolution equations are summarized for several well- 
known models in table [T] Model results are also quantitatively sensitive to implementation details such as cutoffs 
or the choice of light-cone vs. energy splitting. However, these implementation specifics do not seem to lead to 
qualitative changes of the model results. A very instructive review of such effects can be found in IfTTl . 

3. Physics assumptions in jet quenching models 

In addition to implementing different approximations and solution strategies, models also show genuine differ- 
ences with regard to assumptions about the physics underlying parton-medium interaction, and this leads to various 
testable consequences. In particular, the pathlength dependence of the mean energy loss AE of a parton traversing a 
constant medium is parametrically different dependent on what one assumes to be the relevant physics (note that AE 
is not as such a very well-defined concept for an in-medium shower, but it can be extracted by tagging a quark flavour 
as done e.g. in (6)). 

For any incoherent mechanism, the mean energy loss largely tracks the number of scatterings, which in a constant 
medium are proportional to the length, i.e. AE ~ L. For radiative processes, there is a decoherence time involved and 
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hence generic arguments can be made (e.g. [ 1 8 1) that AE ~ L 2 . It is however known that this holds only for sufficiently 
large parent parton energy and finite energy corrections quickly change this to a linear dependence lfl9l l20l I2TI . In 
strong coupling scenarios motivated by AdS/CFT approaches to the medium dynamics, a longitudinal drag rather than 
transverse random momentum transfer lead to induced radiation, which changes the above argument dimensionally, 
leading to AE ~ L 3 [22], to which finite energy correction have so far not been obtained, but would be expected to 
weaken the dependence. 

This discussion is somewhat correlated to the way the medium is implemented in models as discussed in the 
previous section, as models which include explicit interactions with perturbative medium quasiparticles typically find 
50% elastic (incoherent) energy loss for reasonable values of the strong coupling a s ~ 0.3 [ 12, 23]. In contrast, many 
older leading parton radiative energy loss (e.g. I9l [l6l ) start from the assumption of static scattering centers in the 
medium and hence do not have any incoherent component. 

Yet a different source of pathlength dependence is the idea that the minimum virtuality scale down to which a 
shower can evolve in-medium should be constrained by uncertainty arguments to Qq = y/E/L [8 |. This leads to a 
non-linear pathlength as well as additional energy dependence of the mean energy loss [20|. 

It is important to note that in a real fluid-dynamical background the actual pathlength dependence is drastically 
changed by the spacetime evolution of the medium density. In this sense, expressions like 'quadratic pathlength 
dependence' are to be understood as labels what a model would do if applied to a constant medium, not as descriptions 
of what is experimentally measurable. In the following, we aim to establish experimental tests sensitive to these 
different physics assumptions and implementation differences such that the data can be used to discriminate between 
models. 

4. Single hadron suppression 

In general, the Pj dependence of the single inclusive hadron suppression factor Raa is more driven by generic 
pQCD effects such as the functional form of the produced parton momentum spectrum than by jet quenching model 
specific effects (see [28] for a discussion). As a result, within the uncertainty associated with the choice of the bulk 
medium evolution model, models tuned to RHIC data tend to extrapolate reasonably well to LHC conditions (see 
Fig# 

PbPb 2.76 ATeV, 0-5% centrality 



PbPb\'5„ N = 2.76 TeV 




Figure 1 : Pre- and postdictions for the nuclear suppression factor Raa in 2.76 ATeV central Pb-Pb collisions comparison with data from CMS 1 29 ] 
and ALICE (30j in various models I24ll25ll26ll271 . 

The important exception is the strong coupling scenario denoted as AdS which overquenches significantly. This 
can be understood as follows: Due to dimensional reasons, the L 3 pathlength dependence requires a scaling of the 
quenching effect with the medium temperature ~ T A whereas all other scenarios scale with the approximate medium 
density ~ T 3 . Such a strong increase of the quenching power is however not indicated by the data, essentially ruling 
out this particular application of AdS/CFT ideas to jet quenching. 
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Raa is, however, not expected to be a very sensitive probe ll28l — due to the need to average over the full medium 
geometry and to fold with the pQCD parton spectra, only a narrow region of the medium-modified fragmentation 
function (MMFF) is probed by the observable, and thus models predicting a very different MMFF can result in nearly 
the same R AA . One way to constain models better is to study the dependence of R AA on the angle of the hadron with 
the reaction plane <p. In this way, the pathlength dependence of the parton-medium interaction model in combination 
with a particular medium evolution model is tested. Such investigations have established two crucial insights ll28ll3Tll : 

• There is about a factor two uncertainty associated with the choice of the medium model provided that the 
medium model is a fluid-dynamical model constrained by bulk observables. One implication is that R^iip) 
has the capability to discriminate between fluid-dynamical models, i.e. it is a true tomographical measurement. 
However, this result also indicates that any results obtained with a schematical medium evolution model not 
constrained by bulk data (such as the Bjorken cylinder) should be disregarded, as they have easily a factor 10 
systematic uncertainty. 

• Within the factor two uncertainty, most scenarios describe the data reasonably well. The one exception which 
fails by a huge margin (factor 6) is a linear pathlength dependence [32]. Based on comparison with the data, 
any component with linear pathlength dependence must be smaller than about 10% [28, 32 33 1. This does not 
only disfavour scenarios with a 50% elastic energy loss component, but also any radiative scenario with finite 
energy corrections. The only realistic pathlength-dependence generating option viable with the combined data 
is hence the prescription for the minimum in-medium virtuality Qq — y/E/L as suggested in [8 1 (both radiative 
energy loss from an infinite energy parton and an AdS I? dependence would work with the data, but the first is 
not a sufficiently realistic scenario and the second alternative does not agree well with the extrapolation from 
RHIC to LHC as discussed above). 

5. Hard dihadron correlations 

The normalized dihadron away side correlation strength I AA is an interesting observable, since the requirement 
of a hard near side trigger adds a series of biases on parton type and kinematics while the away side fragmentation 
pattern has no additional bias [34]]. As a result, the MMFF is probed in more detail in such measurements. 

This allows to study the redistribution of longitudinal momentum inside a shower in quite some detail. The two 
main physics scenarios outlined above expect rather different behaviour: In a direct, elastic energy transfer into the 
medium, longitudinal momentum lost from the jet is carried by the medium at thermal momentum scales (i.e. not 
in hard modes). In contrast, in a radiative energy loss scenario, longitudinal momentum is carried by additional soft 
gluon radiation, part of which is still 'hard' as compared to a thermal momentum scale. Dihadron correlation are able 
to discriminate between these patterns. 



AuAu 200 AGeV 0-5% centrality 



0.8 



• STAR data 

• YaJEM-D 

• YaJEM-DE 



elastic energy transfer 



ASW 




8 



1.2 



0.4 



0.6 



0.8 



Figure 2: Away side normalized dihadron correlation strength Iaa as compared with leading parton energy loss models and in-medium shower 
models I34II35I . 
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This is illustrated in Fig.|2]where the away side correlation strength as a function of zt — E hadron /E, r i gger for central 
200 AGeV Au-Au collisions and a 8-15 GeV trigger is shown. In leading parton energy loss calculations, energy is 
implicitly assumed to be transferred into the medium and to reappear at thermal momentum scales (below the zt range 
of the measurement). As a result, the curves bend downward at low zt which is not in agreement with the trend seen 
in the data. In contrast, in in-medium shower pictures where the enhanced production of subleading jet fragments is 
explicitly treated, the curves reflect the enhanced soft gluon production by bending upward. 

However, under the assumption that all energy re-appears in soft gluon production, the computation overshoots 
the data. The tend seen in the data can only be reproduced under the assumption that about 10% of the energy is 
directly tranferred into the soft medium [35|. This lower limit for the direct elastic energy transfer agrees nicely with 
the upper limit of the same order as obtained by a study of pathlength dependent observables. Thus, this observable 
clearly shows the limits of leading parton energy loss modelling and also disfavours all scenarios in which there is no 
possibility of direct energy transfer into the medium. 

6. Clustering and suppression of jets 

The clustering of hadron showers into jets by means of a clustering algorithm is designed to suppress the depen- 
dence of observables to soft physics such as hadronization and to provide a good comparison point between hard 
pQCD calculations and measurements without the need to model physics close to Aqcd- As a result, clustering 
suppresses many effects of a medium modifiaction since these take place at a scale T ~ Aqcd- To provide a con- 
crete example, while the medium-induced emission of an almost collinear gluon leads to a modification of the leading 
hadron spectrum since energy has been taken from the leading parton, it does not lead to a modification of the jet spec- 
trum since the emitted gluon is clustered back into the jet. Thus, unless distributions of hadrons inside the clustered 
jets are considered, jets are significantly less sensitive to medium modifications than single hadrons l36l . 

There is, however, a two stage mechanism which is able to suppress jets: While hard partons with pj » T 
are kinematically very robust against interactions with the medium and cannot easily be scattered out of the jet, 
the interaction can nevertheless induce the emission of additional soft gluons at a thermal energy scale. Such soft 
gluons however are not kinematically robust and can be scattered easily out of the jet cone, leading to jet energy loss. 
Since gluons at a thermal scale are indistinguishable from medium gluons, the energy radiated into these modes is 
quickly thermalized and flows in hydrodynamical excitations to large angles. This is thus a very generic mechanism 
independent of the specific in-medium shower physics. The characteristics of the parton-medium interaction are then 
only apparent from semi-soft gluons above the thermal scale. Such a scenario is well in line with jet quenching as 
characterized by RHIC observables discussed before. 

This mechanism of jet suppression has recently been denoted as frequency collimation |37] but has been imple- 
mented into MC and observed earlier in 1 38 1 . Since it relies on extra available phase space for rescattered soft gluons, 
it can not be observed in purely probabilistic approximations of the jet-medium interaction. 

7. Jet observables 

The first LHC observable involving reconstructed jets has been the dijet energy imbalance Aj — f^jr 2 or simply 
the ratio Eji/Etx of reconstructed jet energies Ej\ on the near side and £7-2 on the away side 13911401 . For the reasons 
mentioned above, clustering largely removes the sensitivity to specific characteristics of the jet quenching model, and 
as a result the data can be described in various models making somewhat different assumptions about the dynamics of 
parton-medium interaction and medium geometry ll36l I4T1 l42l l43l . The trigger energy dependence of the imbalance 
is somewhat more constraining, it probes the medium-induced broadening of the jet [44 1 (see Fig. [3]). 

A greater challenge to models is posed by the CMS observation that the longitudinal momentum distribution of 
reconstructed jets, when plotted as a fraction of the reconstructed jet energy, appears unchanged from vacuum even for 
highly imbalanced events B31 . This is a highly unexpected finding for any model which implements the interaction 
with the medium probabilistically and hence uses a modified set of splitting kernels P 1 , (z) to build the fragmen- 
tation function, as in such a model the MMFF would be self-similar, but not show the same functional form as in 
vacuum. It is however a natural outcome for models which break scale invariance only below a fixed momentum scale 
0(T), since in such frameworks the splitting kernels generate a self-similar fragmentation function with the same 
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Figure 3: Dijet asymmetry as measured by ATLAS 1 39 1 and CMS 1 40 1 as compared with medium induced radiative energy loss scenarios 1 4 1 44] 

functional form as in vacuum above the breaking scale. The flattening of Iaa(zt) for zt > 0.5 in Fig.[2]is a manifes- 
tation of the same dynamics. From the RHIC data, significant deviations from the vacuum dynamics are expected to 
occur around 2-3 GeV, below the range of the CMS measurement. Thus, the CMS finding of an apparently unmodified 
fragmentation pattern strongly disfavours an probabilistic implementation of medium effects in the splitting kernel, 
but is not unexpected for other models, indeed almost unmodified jet properties above a momentum cut of 4 GeV with 
significant modifications below had been predicted in [38 1. 

8. Summary 

Our knowledge of the viability of various models, given the available data and aking use of the full uncertainty 
range given by the choice of the hydrodynamical evolution model can be summarized as in Table [2] 





R AA ((f>)@RmC 


R AA @LHC (P T ) 


7 AA @RHIC 


7 AA @LHC 


Aj 


E T 2/E T1 (E, r i g ) 


elastic 


fails 


works 


fails 


fails 


works 


fails 


elMC 


fails 


fails 


fails 


not tested 


not tested 


not tested 


ASW 


works 


fails 


marginal 


works 


N/A 


N/A 


AdS 


works 


fails 


marginal 


works 


N/A 


N/A 


YaJEM 


fails 


fails 


fails 


fails 


works 


works 


YaJEM-D 


works 


works 


marginal 


marginal 


works 


works 


YaJEM-DE 


works 


works 


works 


works 


works 


works 



Table 2: Viability of different parton-medium interaction models tuned to the Pj dependence of Raa in 200 AGeV Au-Au collisions given various 
data sets under the assumption that the best possible hydrodynamical evolution scenario is chosen. The various labels refer to: elastic [33], elMC 
(321 . ASW Q]D, AdS (22), YaJEM (DID, YaJEM-D (25), YaJEM-DE (35). 

It is immediately evident that the combination of all observables cuts the available model space down much more 
than any single observable. Most of the constraints are provided by the combination of the P T dependence of R^ 
at LHC which probes the dynamics of leading partons in combination with a dihadron correlation I AA which is a 
probe of subleading fragment dynamics. Current jet measurements do not have strong constraining power beyond 
this, however they do highlight the need to go beyond a purely probabilistic implementation of medium modification 
and also provide constraints for transverse dynamics of showers. A number of conclusions with regard to the nature 
of jet-medium interaction can be drawn from this: 

• The combined data is compatible with a fairly standard pQCD picture of medium-induced radiation and a 
subleading component of elastic interactions leading to energy transfer into the medium 

• There is no evidence for exotic scenarios, the data neither indicate an L 3 dependence of energy loss as suggested 
by some variants of AdS/CFT inspired models nor a modification of the hadronization stage. 
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• Data from RHIC and LHC consistently indicate that some moderate amount of energy is transferred from the 
hard parton directly into the non-perturbative medium rather than into subleading jet fragments. However, the 
amount of direct energy transfer is constrained by pathlength-dependent observables and found to be signifi- 
cantly smaller than if the medium can be described as a near-ideal quark-gluon gas. This may imply that the 
observed degrees of freedom in the medium are massive or correlated quasiparticles. 

To make these conclusions more quantitative, further high-statistics multi-differential measurements are needed 
from the experimental side, followed by systematic multi-observable studies from theory. In particular, the systematics 
of h-h, jet-h and y-h correlations (with decreasing geometry bias) as a function of the reaction plane angle for various 
system centralities is expected to be a sensitive probe of both the parton-medium interaction mechanism and the 
medium density evolution. Once these questions are reliably addressed, the precise nature of the mechanism by which 
the energy deposited from hard partons is carried by the medium presents itself as the next goal. 
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